Continuous quantum phase transition between an antiferromagnet and a 
valence-bond-solid in two dimensions; evidence for logarithmic corrections to scaling 
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The antiferromagnetic to valence-bond-solid phase transition in the two-dimensional J-Q model 
(an S = 1/2 Heisenberg model with four-spin interactions) is studied using large-scale quantum 
Monte Carlo simulations. The results support a continuous transition of the ground state, in agree- 
ment with the theory of "deconfined" quantum criticality. There are, however, large corrections to 
scaling, of logarithmic or very slowly decaying power-law form, which had not been anticipated. This 
suggests that either the SU(A) symmetric noncompact CP'^"^ field theory for deconfined quantum 
criticality has to be revised, or that the theory for N — 2 (as in the system studied here) differs 
significantly from A — >■ oo (where the field theory is analytically tractable). 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Mg, 75.40.Cx 



Valence-bond solid (VBS) states of two-dimensional 
(2D) quantum spin systems have been studied for more 
than two decades [ij and have recently come into renewed 
focus with the theory of "deconfined" quantum critical- 
ity (DQC) 0,0] I which describes the transition between 
an antiferromagnetic (AF) and a VBS ground state in 
terms of deconfinement of spinous. In addition to the in- 
terest in such AF-VBS transitions in condensed matter 
physics, there are also intriguing connections to decon- 
finement in gauge theories in particle physics [4| . To test 
the validity of the DQC scenario, and to obtain quantita- 
tive results for, e.g., predicted unusual critical exponents, 
unbiased numerical studies of quantum spin hamiltonians 
with AF-VBS transitions are necessary. 

The "J-Q" model was introduced recently Q as an 
SU(2) symmetric spin system realizing the 2D AF-VBS 
transition, following earlier work on related U(l) sym- 
metric models [1, S]. It combines the standard Heisen- 
berg antiferromagnet with four-spin interactions, which 
lead to local correlated bond singlets (valence bonds) 
and reduce the amplitudes of the longer valence bonds 
required in an AF state. The J-Q model is free 
from "sign problems" \^ , which prohibit quantum Monte 
Carlo (QMC) studies of frustrated spin systems such as 
the J1-J2 Heisenberg model on which much of the 
past computational (exact diagonalization) research on 



VBS sates was focused. While series expansions [11 1 
around various candidate states can give some insights, 
QMC methods I4I, when applicable, are the only un- 
biased tools for studying 2D quantum phase transitions 
(in contrast to one dimension, where the density-matrix- 
renormalization-group method [H is apphcable) [3]. 
Being sign problem free, the J-Q model (and generaliza- 
tions of it [15| ) have opened up new avenues for explor- 
ing magnetically quantum-disordered states and quan- 
tum phase transitions. 

In this Letter, a large-scale, high-precision QMC study 
of the AF-VBS transition in the J-Q model is presented 
in order to further test the he DQC theory, and to settle 



discrepancies between previous studies 

The 

main point of contention is the order of the transition. 
In the DQC theory, it was argued^that AF-VBS tran- 
sitions are generically continuous [2] and that the criti- 
cal point for SU(A'^) spins corresponds to a non-compact 
(NC) CP^-i field theory f3]. This is at odds with the 
long-standing Landau-Ginzburg paradigm, where a di- 
rect transition between two states breaking unrelated 
symmetries should be first-order (except at fine-tuned 
multi-critical points) . Ground state [1, [11] and finite- 
temperature 'l6i| QMC studies of the J-Q model show 
scaling behavior in good agreement with the DQC the- 
ory, including a dynamic exponent z = 1, a rather large 
anomalous dimension Jospin ~ 0.35, and an emergent U(l) 
symmetry in the VBS phase (which in the theory is asso- 
ciated with spinon deconfinement). On the other hand, a 
QMC finite-size analysis by Jiang et al. would, if correct, 
require a first-order transition [1^. A weakly first-order 
AF-VBS scenario has been elaborated by Kuklov et al. 



19|, based on results for a lattice model claimed to 
realize the NCCP^ action, but other studies of the action 
have reached different conclusions 

Here it will be shown that the claimed first-order sig- 
nals in the study by Jiang et al. [17] can be attributed 
to over-interpretations of QMC data affected by signif- 
icant systematical and statistical errors. The results to 
be presented below were obtained with the stochastic se- 
ries expansion (SSE) method 0,[2l|, which is a finite- 
temperature QMC method free from systematical errors. 
There are no indications of a first-order transition, even 
in systems of space-time volume 20 times larger than in 
[13]. However, the data are now of high enough quality 
to detect logarithmically weak deviations from the scal- 
ing forms expected at a z = 1 critical point. Logarithmic 
corrections are well known consequences of marginal op- 
erators at criticality, which, although they have not been 
predicted theoretically in this case (in large- TV treatments 
of the NCCP^^i theories @,|2^,|23|), cannot a priori be 
ruled out for iV = 2. A first-order transition would lead 
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FIG. 1: (Color online) Binder cumulant for the sublattice 
magnetization as function of the coupling ratio for different 
system sizes at inverse temperature (3 — L. 

to much more dramatic deviations from z — 1. 

Turning now to a quantitative discussion of the calcu- 
lations, the J-Q hamiltonian '3] can be written as 



Q J2 

(ijki) 



ij Ckl , 



(1) 



where Cy is a bond-singlet projector for S* = 1/2 spins; 
Cij = 1/4 — Si • Sj . In the J term ij are nearest neighbors 
on the square lattice, while ij and kl in the Q term are on 
opposite edges of a 2 x 2 plaquette. Lattices oi N = L? 
spins with periodic boundaries are used. Assuming z = 1 
(based on previous work the inverse temperature 

P = Q /T is taken proportional to L for finite-size scaling; 
(3 = L and /3 = i/4 will be considered for L up to 256. 
Calculations for T/Q > 0.035 are also carried out for 
systems sufficiently large, up to L = 512, to give results 
in the thermodynamic limit. 

The focus here will be on magnetic properties. The 
staggered magnetization is computed along the z 
(quantization) axis. To extract the critical coupling ratio 
{J/Q)c, and to address the issue of a possible first-order 
transition, consider first the Binder cumulant [24], 



U2 = - 1-- 



1 



3 (toL)^ 



(2) 



which is defined so that U2 and C/2 — 1 in an AF 
disordered and ordered state, respectively, when L ^ 00 
(stemming from a Gaussian distribution of around 
|ms| = and a (5-function at \rfis\ > 0, respectively). The 
factors in ([2]) correspond to rn^z being one component of 
a three-dimensional vector fris. At a continuous transi- 
tion, curves oi U2 versus J/Q for different system sizes 
should intersect at the critical coupling, where normally 
< t/2 < 1 [13] ■ At a first-order transition, on the other 
hand, U2 ^ —00 when L ^ 00 
tribution with peaks at both \rris 



2J] , following from a dis- 
> and \ms \ = when 
the ordered and disordered phases coexist (with weight 
transferring rapidly between the peaks as the transition 
is crossed for large finite V). It should be noted that U2 
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FIG. 2: (Color online) Critical couplings extracted from the 
crossing of Ui(IS) and U2(L/2) and from the winding number 
criterion Po = 1/2 in systems with P = L and L/A. 



can be negative also at a continuous transition [2J, |25| — 
only a divergence signals a first-order transition. 

As seen in Fig. [1] in the J-Q model there are no signs 
of C/2 becoming negative. The curves intersect at a point 
which moves very slowly toward larger J/Q with increas- 
ing system size. The critical coupling for L ^ 00 can be 
extracted by extrapolating the crossing points for sys- 
tems of size L and L/2, as shown in Fig. [2j 

Fig. [2] also shows results for the size-dependent critical 



coupling su gges ted by Kuklov et al. [18| and used by 



Jiang et al. [17|. It is based on the winding numbers. 
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(3) 



where Ja{p), a = x^y, is the spin current in lattice direc- 
tion a at location p in an SSE configuration containing 
n operators [2^. In the case of the J-Q model, these 
currents take the values Jaip) G {0,±1,±2}. The "tem- 
poral" winding number is essentially the magnetization; 



Wr = 2M^ 



N 

E 



Sf. 



(4) 



The squared winding numbers are related to two impor- 
tant thermodynamic quantities; the spin stiffness. 



and the uniform magnetic susceptibility. 
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(5) 



(6) 



For L ^ 00 and T — >• 0, in a magnetically disordered 
(here VBS) phase ps ^ Q and x — 0, while in the AF 
phase ps > and x > 0. A possible definition of the 
transition point (for finite L and /3) is the coupling at 
which the probability Po of all the winding numbers be- 
ing zero is 1/2 (or any fixed fraction) [l3|. Fig. [2] shows 
results obtained by interpolating Pq for several J/ Q val- 
ues. They extrapolate to the same {J/Q)c ~ 0.0445 as 
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FIG. 3: (Color online) The uniform susceptibility divided by 
the temperature in the neighborhood of the critical point. 
The solid curve is of the form x/T = a + 61n(Q/r) + cT^ . 

the Binder cumulant crossings, but the corrections are 
larger. Note also that the size dependence varies sig- 
nificantly with the aspect ratio P/L. The results do not 



agree well with those of Jiang et al. [17[ , although ps and 
X agree reasonably well for the system sizes available for 
comparisons. It is possible that Pq is more sensitive to 
the Trotter approximation used in |T3|- Note also the 
non-monotonic size dependence in Fig. [2j An \n{L)/L^ 
convergence of {J/Q)c was cited in [l^l as a sign of a first- 
order transition. The data fits were, however, based on 
only three system sizes. The behavior for larger lattices 
is clearly different. 

The L ^ oo critical value {J/Q)c ~ 0.044 is marginally 
higher than in previous studies. In particular, fitting the 
expected z = 1 form x ~ P of the susceptibility at T > 
(L ^ oo), Melko and Kaul found (J/Q)c « 0.038. At 
higher J/Q they found x = a + bT, as expected in the AF 
phase. Fig. [3] shows x/T down to temperatures less than 
half of the lowest T considered in At J/Q = 0.04, 
while x/T is roughly T-independent for 0.05 < T/Q < 
0.2, there is a drop at lower T, consistent with a spin- 
gapped phase. Close to the critical point there is no pure 
X oc T dependence at low T; instead the data exhibit a 
slow divergence, x/T w a -f bhi(Q/T). The fanning-out 
of the data suggests that the logarithmic form is a critical 
separatrix between the expected T ^ behaviors in the 
VBS and AF phases. 

Another indication of logarithmic corrections comes 
from the total squared winding number, 



{W^) = 2Pp, + ^x, (7) 



for which Jiang et al. claimed an asymptotic linear diver- 
gence at the transition as would be expected when 
AF and VBS phases coexist at a first-order transition. 
Fig. m shows the results of the present study. While (W^) 
indeed grows with L, it does so very slowly, consistent 
with a logarithmic divergence. There is no plateau fol- 
lowed by a linear divergence — that conclusion [l7j seems 
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FIG. 4: (Color online) Size dependence of the total winding 
number at the Po{L) = 1/2 point for (3/L and 1/4. The 
data are shown on log-lin (main panels) and lin-lin scales (in- 
sets). The solid and dashed curves are fits to forms a-|-61n(L) 
and c — ciZ/~°'^, respectively. 



to be based on an over-interpretation of noisy data. 

In principle, it is not possible to distinguish between 
a logarithm and a conventional scaling correction ^ L^'^ 
with a very small > 0. Fig. 2] shows fits with w = 0.1 
along with the logarithmic form — when a; — > the two 
forms coincide exactly. This comparison shows that if 
the corrections are conventional, then lj < 0.1. This is 
true also for the uniform susceptibility (Fig. [3]). 

Consider now the stiffness [not combined with x as in 
([7])]. At a conventional z = 1 critical point ps ^ 1/L. 
In the present case the drift in crossing points of psL 
curves for different L is larger than what is normally 
(26| expected, but can be compensated by a logarithm, 
PsL/\n{L/LQ), as shown in Fig. [5] For L > 48 the curves 
intersect at a point, giving {J/Q)c = 0.0447 ± 0.002, 
in agreement with all the other results discussed above. 
Scaling fits away from the critical point give a correlation 
length exponent u « 0.6, but this is without considering 
possible corrections also to the conventional L^^'^ scaling. 
It is difficult to include logarithmic corrections in quanti- 
ties where the leading exponent is not known, in contrast 
to Ps and X where z = 1 governs the leading behavior. 

The conclusion of this study is that the AF-VBS tran- 
sition in the J-Q model is continuous, but with signifi- 
cant corrections to the z = 1 scaling that have not been 
discussed previously. The corrections appear to be loga- 
rithmic, although conventional scaling corrections ~ L~'^ 
with w < 0.1 cannot be ruled out based on the numerical 
data alone. Regarding the possibility of a very weakly 
first-order transition, it should be noted that rigorous 
proofs of continuous phase transitions are only available 
for a small number of exactly solvable models, yet accu- 
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FIG. 5: (Color online) Scaling of the spin stiffness with a log- 
correction (I/O = 0.9) for P = L systems. The curves show a fit 
to a common polynomial /[(J — Jc)!/^'^"] with Jc/Q = 0.0447 
and u — 0.59 (not including the L = 32 data). 

mulated numerical evidence of scaling (and experiments 
on natural systems), along with non-rigorous analytical 
calculations, have established a consensus that critical 
points are ubiquitous. The system volumes /3L^ used 
here for the J-Q model are similar to those in contem- 
porary classical Monte Carlo simulations [27]. In the ab- 
sence of any concrete signals of first-order behavior, the 
transition must therefore be regarded as continuous. 

The scaling corrections will hopefully stimulate further 
field-theoretical work to explain them. Scaling anomalies 
that could be logarithmic have been seen in Monte Carlo 
studies of the NCCP^ action but it has also been 
claimed that this action always leads to a first-order tran- 
sition [l9[ (in which case a different field-theory for the 
J-Q model would have to be found). Marginal operators 
leading to logarithms appear in systems at their upper 
critical dimension, but this is not applicable here. Loga- 
rithmic corrections have been previously found in gauge 
field theories with fermions [28]. On the other hand, con- 
ventional power-law corrections due to irrelevant opera- 
tors are always expected, but here the subleading expo- 
nent uj would have to be very small, which has not been 
anticipated (although the dangerously irrelevant opera- 



tor causing the VBS has a small scaling dimension [15 1 
and is a potential source of a small lj). Studies of the 
SU(A^) generalization of the J-Q model would be useful 
to determine whether iV = 2 is a special case. QMC cal- 
culations have already been carried out for = 3 and 4 
15{ , but the quantities discussed here have not yet been 
investigated. 

A consequence of the findings presented here is that the 
anomalous VBS transition in U(l) symmetric systems Q 
should be re-evaluated. Scaling deviations very similar to 
(but stronger than) those in the J-Q model were found, 
which in jisl [l9| was interpreted as a first-order transi- 
tion. Considering scaling corrections, this class of models 
as well may in the end have continuous transitions [6] . 
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